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Abstract. The Skyrme-Hartree-Fock-Bogoliubovcode HFBTHO using the axial (2D) Transformed 
Harmonic Oscillator basis is tested against the HFODD (3D Cartesian HO basis) and HFBRAD 
(radial coordinate) codes. Results of large-scale ground-state calculations are presented for the SLy4 
and SkP interactions. 



1. INTRODUCTION 

The code HFBTHO |l| solves the self-consistent HFB equations by using the axial 
(2D) Transformed Harmonic Oscillator (THO) basis 10], which allows for a correct 
treatment of the single-quasiparticle wave function asymptotics. As discussed recently 
the THO technique is a method of choice for performing massive nuclear structure 
calculations including weakly bound systems. In order to fully test the formalism, in the 
present study we present results obtained with the axial (2D) HFBTHO (vl.64) code 
compared to those obtained with two other codes: HFODD (v2.08i) fl], which uses a 
Cartesian (3D) Harmonic Oscillator (HO) basis and spherical (ID) HFBRAD which 
uses a lattice of points in the radial coordinate. 

In Ref. im we have published the first complete mass table of even-even nuclei 
obtained by using the THO method for the SLy4 Skyrme force Here we discuss one 
specific improvement of the method, and also present new results obtained with the SkP 
Skyrme force 10] ■ More details, including downloadable tables of ground-state proper- 
ties can be found atjhttp : / / www . f uw . edu . pl/~dobaczew/ thodri/thodri . html 



2. TESTS 

In this section, we discuss results of two numerical tests. First, by switching off the Local 
Scaling Transformation (LST) of THO, we run HFBTHO in the axial HO basis and test 
it against HFODD. For a given Skyrme interaction and zero-range, density-dependent 
pairing force, both codes should give exactly the same results. Since technical details of 
the inner structure of both codes are completely different, such calculations constitute 
an extremely stringent test of both codes. 



Second, by switching the LST on, we could test the code HFBTHO against the 
spherical code HFBRAD fsj] . Here, results of both codes cannot be exactly identical, 
because the phase spaces in which the solutions are obtained are significantly different. 

Table [H displays the results of test calculations performed for the SLy4 Skyrme in- 
teraction [6] and for the mixed zero-range pairing force fsl]: y(r)=Vb(l ^p(^)/po) for 
po=0.32fm^-^. The cutoff energy of £t„f=60MeV was used for summing up contribu- 
tions of the HFB quasiparticle states to density matrices For a given phase space, 
the strength of the pairing force Vq was adjusted so as to reproduce the experimen- 
tal neutron pairing gap in ^^^Sn. The resulting values are Vo=— 285.88, —284.10, and 
-284.36 MeVfm^ for the HO (THO) bases of 680 and 3276 states, and for the ra- 
dial box of Rbox=^^ fm, respectively. The radial HFBRAD calculations were performed 
with 300 points (i.e., the Ar=0.1fm grid spacing), and the wave functions were in- 
cluded up to j,„ax=^9/2. We checked that even with j„„a=33/2, all energies were sta- 
ble within 1 eV. The nucleon-mass and elementary-charge parameters were fixed at 
^V2m=20.73553MeVfm2 and 1.43 997 g MeVfm, respectively. 

Table [T] displays the following quantities: A^o is the maximum number of the HO 
oscillator quanta included in the basis (for the deformed basis we give the numbers 
of quanta in the perpendicular (A''j^) and axial (A^,) directions); Nst is the number of 
the lowest deformed HO states included in the basis; Nn'^ and Np'^ are the numbers 
of (doubly degenerate) neutron and proton quasiparticle states with equivalent single- 
particle energies |2] below the cutoff energy Ecut', ^± and are the oscillator constants 
in the perpendicular and axial directions; A„ and Xp are the neutron and proton Fermi 
energies, which, for vanishing pairing correlations, are taken as the s.p. energies of the 
last occupied states; A„ and Ap are the average pairing gaps O]; i?„ and Rp are the 
rms radii; Q„ and Qp axe the quadrupole moments (2z^ —x^ —y^) \ £n^ and ef,^ are the 
s.p. energies of the most bound neutron and proton states; Lf^ and are sums of the 

canonical energies weighted by the corresponding occupation probabilities; En'^"'^ and 

Ep""^ are the pairing energies; E^"'' and Ep^" are the kinetic energies; Ecen and Eso are the 
energies corresponding to the central and spin-orbit parts of the Skyrme energy density 
functional; E^ir and Eexc are the direct and exchange parts of the Coulomb energy; and 
Estah is the stability energy characterizing the level of self-consistency. In the code 
HFODD, 

Estah is estimated from the sum of s.p. energies ^; in the code HFBTHO 
Estab is estimated from the maximum difference of all matrix elements of s.p. potentials 
calculated in two consecutive iterations; and in the code HFBRAD it is calculated as a 
variance of the total binding energy, Ejot, over the last five iterations. 

Calculations for ^^^Pb yield a spherical solution with vanishing pairing gaps. HF- 
BTHO and HFODD give the total binding energies that differ by 627 eV, and this dif- 
ference can be (primarily) traced back to the direct Coulomb energy. We have checked 
that without the Coulomb interaction, this difference decreases to 202 eV. The axial- 
basis HFBTHO calculation gives a very small total quadrupole moment of 39 /ib. This 
suggests that the THO basis generates a slight deviation from the spherical symmetry 
due to a different numerical treatment of z- and ^-direction. In this respect, HFODD 
calculations should be considered more accurate. 

Calculations for ^^^Er performed within a spherical HO basis, b±=bv, yield a well- 
deformed and weakly paired prolate ground state. Here, the total binding energies and 



TABLE 1. (Color online) Benchmark results of the HFB calculations performed for the SLy4 interaction and mixed 5 pairing. All energies are in MeV, 
lengths in fm, and quadrupole moments in barns. Boldface colored digits differ between the HFBTHO and HFODD/HFBRAD calculations. See text for 
details. 
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FIGURE 1. (Color online) Left: differences in Efot obtained in HFBTHO by using LST based on 
RHO or LAM conditions. Right: differences between Eioi obtained in THO and HO bases. Calculations 
were performed using the SLy4 interaction with volume pairing and 20 oscillator shells. Lipkin-Nogami 
method followed by the exact particle-number projection was used to correct for the particle number 
nonconservation in HFB. 



quadrupole moments obtained within HFBTHO and HFODD differ only by 32 eV and 
5 /ib, respectively. When the same calculation is performed in a deformed HO basis, 
b±^bz, the differences grow to 5.1 keV and 207 /ib, respectively. Again, without the 
Coulomb interaction, the difference in the total binding energy is only 96 eV. It is seen 
that by employing the deformed basis, the binding energy decreases, as expected. 

Comparison with the coordinate- space code HFBRAD for ^^^Sn shows that Etot in 
HFBTHO is correct up to 14keV for No=25. However, the kinetic energy still differs 
by as much as 221 keV, which is compensated by a similar difference in the interaction 
energy. Within the HO basis and A^o=25, the corresponding differences are larger: 41 and 
337 keV. The analogous differences obtained for A^o=20 are 142 and 1103keV (THO), 
and 152 and 964 keV (HO), respectively. Nevertheless, the above comparison shows that 
the No=20 calculations yield Etot with a precision of a couple of hundred keV. 



3. MASS TABLES 

The LST employed in Ref. f3\ was based on HO densities corrected in the asymptotic 
region by the contribution from the lowest-energy quasiparticle. Since a common LST 
has to be carried out for both neutrons and protons, for each nucleus one is forced to 
make a decision whether the LST is to be based on neutron or proton density. In Ref. 10] 
we used a prescription (referred to as LAM) that the neutron densities were used for 
A„>Ap and vice versa. In this work, we use the condition pn{Rmin)>Pp{Rinin), where 
Rmin is the point where the neutron or proton logarithmic density has a minimum as a 
function of r. In practice, the above condition, dubbed RHO, does not depend on whether 
the neutron or proton 7?„„„ is considered. 

In Fig. [U (left panel) we show the differences in Etot obtained in HFBTHO by using 
the LST condition employing the Fermi energies (LAM) [3^ or the densities (RHO). One 



can see that in the majority of neutron-rich nuclei both prescriptions lead to identical 
results. However, in many proton-rich nuclei the new prescription decreases binding 
up to about 500 keV, and for some medium-mass proton-rich nuclei the RHO method 
decreases binding by up to lOOkeV. This latter effect is due to a better description of 
asymptotics in the pairing channel, which leads to extended pairing fields and reduced 
pairing energies llioll . The right panel of Fig. [l] shows differences in Etot obtained in 
THO and HO bases. In most nuclei, by using the THO basis, one obtains a small 
energy gain of up to lOkeV. This grows to ~500keV for the very neutron-rich systems. 
Again, in lighter nuclei, a better asymptotics may lead to a reduced binding. In fact, 
our results show that improvements in density profiles at large distances cannot be 
treated variationally. First, Etot is quite insensitive to the precise description of nucleonic 
densities in outer nuclear regions. Second, due to the pairing-space cutoff, the pairing 
energy is not reacting variationally on the improvement of the wave function. 

Figures El and El present HFBTHO results obtained with the SLy4 and SkP Skyrme 
forces. It is obvious that without further improvements these traditional Skyrme forces 
describe nuclear masses rather poorly. The rms deviations between calculated and mea- 
sured masses are as large as 3.14MeV for SkP and S.lOMeV for SLy4, respectively, as 
compared to about 0.70 MeV deviations obtained for forces fitted specifically to masses 
(see Ref. lITlll for a review). Moreover, pronounced kinks obtained at magic numbers 
suggest that the quality of the description of (semi)magic and open-shell systems is not 
the same. This may point to a need to systematically include dynamical zero-point cor- 
rections l[l2ll . Work in this direction is in progress. 
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FIGURE 2. (Color online) Ground-state deformations j3 (left) and two-neutron separation energies S2n 
(right) obtained within HFBTHO using SkP (top) and SLy4 (bottom) interactions. 
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FIGURE 3. (Color online) Deviations of ground-state HFBTHO energies from experiment 11311 for 
SkP (left) and SLy4 (right) interactions. Positive values correspond to underbound nuclei. No corrections 
beyond mean field were included. 
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